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Using density matrix equations of motion and a tight -binding band calculation, we predict all- 
optical switching between four metastable magnetic states of (III,Mn)As ferromagnets. This switch- 
ing is initiated non-thermally within lOOfs, during nonlinear coherent photoexcitation. For a single 
optical pulse, magnetization reversal is completed after ~100 ps and controlled by the coherent fem- 
tosecond photoexcitation. Our predicted switching comes from magnetic nonlinearities triggered by 
a femtosecond magnetization tilt that is sensitive to un-adiabatic light-induced spin interactions. 

PACS numbers: 78.47.J-, 78.20.Ls, 42.65.Re 



The goal of THz magnetic switches underlies the en- 
tire field of spin-electronics and challenges our under- 
standing of fundamental non-equilibrium spin processes. 
The reading and writing of bits rely on reversing the 
magnetization direction between "up" and "down" . In 
conventional switching, the magnetization moves out of 
equilibrium via laser heating. A magnetic field then ex- 
erts a torque that reverses the magnetization within few 
ns [HQ. This speed can be improved by using coher- 
ent spin rotation, via precession of the entire memory 
cell around a magnetic field pulse- precessional/ballistic 
switching This pulsed field must have duration of 

at least half the precession period (100's of ps), which 
sets a fundamental limit of the magnetization reversal 
time. This speed is further limited by randomness Q 
or weak precession damping allowing back-switching of 
magnetic elements ( "ringing" ) 0, Q . Faster switching, 
within lOOfs, could be explored by using laser pulses to 
inject spin-polarized carriers [f|. This is important for 
meeting the demand for improved read/write speeds, bit 
density, and reliability of current magnetic devices. 

Magnetic properties of (III,Mn)V ferromagentic semi- 
conductors exhibit sensitive response to carrier density 
tuning via light, electrical gates, or spin currents. This 
holds promise for high-speed magnetic switches that com- 
bine information processing and storage on a single chip 
device with low power consumption [6|. The femtosec- 
ond photoexcitation of GaMnAs revealed distinct tran- 
sient magneto-optical responses: (i) ultrafast decrease of 
the magnetization amplitude [?] within ~ lOOfs (demag- 
netization) [8l 4lOI|. f ii) enhancement of magnetic order 
on ps timescale [ll|, (iii) magnetization re-orientation 
within ~100fs, followed by a distinct ps regime of coher- 
ent precession 1^, 13 1. Such non-equilibrium magnetic 
effects appear to be universal @, The pioneering 

work of Bigot and collaborators J], who observed de- 
magnetization on a ^lOOfs timescale much shorter than 
the spin-phonon relaxation time, thus evolved into a new 



field of femto-magnetism. However, the many-body the- 
ory of femto-magnetism remains controversial [5|, |15[ and 
must ultimately engage the elements of transient coher- 
ence, correlation, and nonlinearity on an equal footing. 
Here we present such a mean field theory and propose a 
nonthermal mechanism for achieving ultrafast all-optical 
magnetization switching in ferromagnetic (Ga,Mn)As. 

We use a single lOOfs optical pulse at ^3eV, in reso- 
nance with the strong peak in the density of states for 
A3— ^Ai interband transitions along the eight equivalent 
1 111} directions of the GaAs Brillouin zone (the A-edge) 
16]. Magnetization switching is initiated within ^lOOfs 
via non-thermal spin manipulation. To describe this, 
we consider the Hamiltonian H (t)=Hb+ H exc h(t)+ H h{t) , 
where Hb=Ho + Hso + Hpd is the standard Hamiltonian 
describing un-excited (Ga,Mn)As bands @, H de- 
scribes the states in the presence of the lattice potential, 
while Hso is the spin-orbit interaction [l7j]. We con- 
sider hole densities ^10 20 cm~ 3 where the virtual crys- 
tal approximation applies @. H pc [ is the mean-field in- 
teraction 0, [l8| of the hole spin with the ground state 
Mn spin So, whose strength /3=24 meV nm 3 was ex- 
tracted from experiment [121 ]. H p d lifts the degeneracy 
of the GaAs bands [ill by the magnetic exchange en- 
ergy A p( j=88meV. Below we describe microscopically a 
photo-induced magnetic anisotropy characterized by the 
energy ratio Eso/^-pd (i?so~350meV is the spin-orbit 
energy), which is not due to thermal 14| effects. For 
this we start by diagonalizing Hb using the Slater-Koster 
sp 3 s* tight-binding Hamiltonian [I?} • We thus obtain 
a basis of valence hole (conduction electron) states cre- 
ated by h^ n (e^ ), where k is the crystal momentum and 



n labels the 20 bands [17[. This static bandstructure 
is modified by the interaction H exc h{t) of the hole spin 
with the photoexcited deviation AS(t) of the collective 
Mn spin S from So 13j, LL8| . Here we treat H exc h (t) at 
the Hartree-Fock level, which conserves the magnetiza- 
tion amplitude, and note that demagnetization |l0( arises 
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from correlations considered in Ref. [19j. Hi,(t) describes 
the interband dipole coupling of the optical field E(t), 
which propagates along [001] and is linearly polarized at 
an angle of ^2° from [100]. Hl leads to nonlinear cou- 
plings between the bands, characterized by the Rabi en- 
ergies d nm k(t)=d nm k cxp[-iuj p t - t 2 /Tp], where w p =3eV, 
r p =100fs, and d nm k=iJ, nm \ s E. The matrix elements /x nm k 
were obtained from the tight-binding parameters of Hb 
as in Ref. [2(|. We derived coupled density matrix equa- 
tions of motion for S, all nonthermal populations, and 
the photoinduced transient coherences between all con- 
duction and valence bands [IH [HI- The dynamics of 
these coherences, neglected within the semiclassical rate 
equation treatment of spin photoexcitation, is important 
for determining the photoexcited hole spin S/j(t). For de- 
tails, see Refs. IH 17, 18, 20| and a future publication. 
The Mn spin is obtained from the equation of motion 



d t S = -/3S x Bh - 7 S x H th 



x d t S, (1) 



where 7 is the gyromagnetic ratio and a the Gilbert 
damping coefficient 2l|. Our mean-field approximation 



misses demagnetization [Ioj |. which can be included phe- 
nomenologically [lij or microscopically [lj|- Here, the 
photoexcited hole spin creates an effective magnetic field 
pulse, r yH p f l (t)=(3sh(t), which transiently modifies the 
magnetic anisotropy and triggers nonthermal switching, 
is obtained microscopically after noting that 



y Skrm^kn^kn) + y S knn' (^kn^kn' ) , 

nk n^n'k 



(2) 

where V is the volume. s£ nm are the hole spin matrix 
elements with respect to the eigenstates of Hi,. They de- 
scribe spin mixing due to H$o and H p d- The first term in 
Eq.@ is roughly proportional to the photoexcited hole 
density. For w p =3eV, it comes from the transient popu- 
lation of high energy states close to the A point; these do 
not contribute to the ground state magnetic anisotropy. 
The second term in Eq.([3]) comes from the coherences 
(h^hkn'), n^n' , photoinduced via Raman processes. 
Here we solved the nonlinear equations of motion that 
give all (/i kn /ikm} at the eight k's with maximum density 
of states (A-point). The mean field arising from all k 
will be considered elsewhere. The dynamics depends on 
H exc h(t) and the nonlinear polarizations (fe-k»ekm)i ob- 
tained from their equations of motion 13 , 18[ . We thus 
describe the effect on Sh(t) of (i) static bandstructure and 
the competition of spin-orbit and magnetic exchange in- 
teractions, (ii) photoinduced interactions that transiently 
change the bandstructure, and (iii) nonlinear coherence. 
We considered dephasing/relaxation times ~30fs [f|. 

The r -point thermal hole Fermi sea Q is clearly dis- 
tinguished from the carriers that create s^(t), which are 
photoexcited at ~3eV. Our main assumption is that, 
within time intervals ^lOOfs, the Fermi sea adjusts adi- 
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FIG. 1: (color online) (a)-(c): S(t) for #=100mT, 
a=0.03 [U, and three pump fluences. _K" c =0.0144meV, 
i('i 1 =0.0025meV, and K UZ =5K C give anisotropy fields 1.7mT, 
0.29mT, and 8.3 mT @|. (d): Angle cos9=S ■ S(t), 
0<O(t)<7r. (e): In-plane component of the non-thermal field 
jHp h =l3(sh) ± ■ (f): Thermal field component H^,. 



abatically to the transient changes in the Mn spin or- 
der parameter and can thus be described in terms of its 
total energy E^[S(i)]. Eh then depends on the instan- 
taneous Mn spin (unit vector S) and the second term 
in Eq.([T]) describes precession around the thermal mag- 
netic anisotropy field jH t h=—dEh/dS [EH]- Symmetry 
22, [23| dictates that, independent of issues surrounding 
the bands close to the bandgap of (Ga,Mn)As (il. Hoj. 



Eh = K c (S x S y +S x S z +S y S z )+K uz S z —K u S x Sy—"fHS z 

(3) 

K c is the cubic 
is the uniaxial constant due 



23] 



as observed experimentally 
anisotropy constant, K uz 
to strain and shape anisotropies, K u describes an in- 
plane anisotropy that may be due to materials issues [23| , 
and H is an external magnetic field along [001]. We ex- 
tract these parameters from experiment [231 ] and neglect 
any transient temperature effects, which enhance our 
predicted anisotropy and give demagnetization 0,0. 
We thus obtain a biaxial magnetic aniso trop y with jour 
metastable magnetic ground states [12j, [24[. Since the 
easy axes are close to y=0 or x=0, wc label these four 
magnetic ground states by AT* and Y ± [I2I ]. 

FiguresQJa), (b) and (c) show S(t) for initial condition 
along A + and three pump fluences. They clearly demon- 
strate magnetic switching controlled nonthermally by the 
pump pulse. Fig. HJa) was obtained for £J=2xl0 5 V/cm, 
as in the experiment of Ref.[l2} (fluence ^7/iJ/cm 2 ). In 
this case, AS is small and excites magnon oscilla- 
tions. These results agree with Ref. [12]. With increasing 
pump intensity, the magnetic nonlinearities of Eq. ([3]) 
kick in. For S=5xl0 5 V/cm, we obtain A— >Y switch- 
ing after ~400ps (to Y^) 1 while for intermediate times, 
the magnetization reverses direction (to X~, FigfT](b)). 
For E=6x 10 5 V/cm, the magnetization reversal is com- 
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FIG. 2: (color online) Nonlinear magnetization dynamics 
(a=0.3, £=6xl0 5 V/cm). (a): Comparison of Mn spin re- 
orientation angles 6(t) and thermal anisotropy field H% h along 
[010] (inset) during ~100ps and ~10ps timescales. (b): Com- 
parison of full calculation (solid line) with result obtained by 
neglecting H exc h(t) (dashed line) and within semiclassical adi- 
abatic approximation of spin photoexcitation (dotted line). 



plete after ~400ps (to X~ , Fig[TJ (c)). FigQJd) shows 
this switching more clearly, by plotting the magnetiza- 
tion reorientation angle 0(f) = cos _1 (S • S) (O<0<7r). 
Our scheme could thus potentially be used to write mul- 
tiple bits in a massively parallel memory with THz speed, 
since one normally reads bits long after writing them. 

Our switching is initiated by the photoinduced non- 
thermal (Hph) field (FigQTe)) and completed by the ther- 
mal (U th ) field (Fig. QJf)). H ph only lasts for ~100fs 
and is enhanced by the external magnetic field H. For 
given intensity, H p h is much stronger for w p =3eV than 
for w p =1.5eV due to the difference in the density of states 
and spin-orbit interaction (l3| . In contrast, the thermal 
field Uth develops as Eh changes due to the Mn spin 
AS(t) photoinduced by s/j. It is much weaker than H p / t 
during the first 100's of fs but dominates over ~100ps. 

We now turn to the origin of the photo-induced nonlin- 
ear behavior demonstrated by FigJlJ By expanding Uth 
to first order in AS, we obtain the linear spin dynamics 
shown by the dotted lines in FigJ2] (a). The difference 
between the linear and nonlinear Mn spin re-orientation 
angle 0(t) is small during fs timescales but grows over 
many ps. To interpret this, the inset of Fig. [5] shows that 
H% h , the component perpendicular to So, grows over 10's 
of ps due to the cubic anisotropy nonlinearity ocSy, as 
S v increases due to precession around H* h . In turn, S z 
(and thus H* h ) increases due to precession around H^ h . 
Switching is triggered by this nonlinear dependence, as 
H t h builds over 10's of ps until the magnetization over- 
comes the energy barrier between the magnetic states. 

FiglUJb)) shows that the time-dependent changes in 
the bandstructure H^, due to the photoinduced inter- 
action H exc h(t), affect sh and 0(t) (compare solid and 
dashed lines). The dotted line in FigJ^fb) also shows 
that the semiclassical Fermi's Golden rule hole spin gen- 
eration rate, obtained by solving our equations within the 
adiabatic approximation while treating the interactions 
perturbatively [25], can also give strong deviations. We 
conclude that our treatment of the coherent nonlinear 
photoexcitation of the interacting system is important 



for modelling the magnetization dynamics. 

In summary, we predict spin switchings triggered by 
coherent nonlinear photoexcitation of (III,Mn)V ferro- 
magnets by a single femtosecond optical pulse. We 
obtained all-optical magnetization reversal and demon- 
strated a four-state magnetic switching functionality re- 
sulting from the interplay between ultrafast coherence, 
nonlinearity, and competing spin interactions. Our pro- 
posed effect should be confirmed experimentally with 
femtosecond magneto-optical spectroscopy and points 
to future possibilities for further reducing the switching 
times by using multiple coherent optical pulses. 
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